New results for 5-point functions 

J. Gluza 1 and T. Riemann 2 * 

1 - Institute of Physics, Univ. of Silesia, Universytecka 4, 40007 Katowice, Poland 

2 - Deutsches Elektronen-Synchrotron DESY 
Platanenallee 6, D-15738 Zeuthen, Germany 



Bhabha scattering is one of the processes at the ILC where high precision data will 
be expected. The complete NNLO corrections include radiative loop corrections, with 
contributions from Feynman diagrams with five external legs. We take these diagrams 
as an example and discuss several features of the evaluation of pentagon diagrams. 
The tensor functions are usually reduced to simpler scalar functions. Here we study, as 
an alternative, the application of Mellin-Barnes representations to 5-point functions. 
There is no evidence for an improved numerical evaluation of their finite, physical 
parts. However, the approach gives interesting insights into the treatment of the IR- 
singularities. 



1 Introduction 

Bhabha scattering, 

e+ + e~ -> e+ + e~ , (1) 

is one of the most important reactions at e + e~ collidersQ At ILC energies, small angle 
Bhabha scattering is dominated by pure photonic contributions and is foreseen as a lu- 
minosity monitor, and large angle Bhabha scattering is also one of the reactions with an 
expected jigh event statistics and with a very clean theoretical Standard Model prediction. 
For these reasons, a NNLO (next-to-ncxt-to leading order) prediction of the complete QED 
contributions and a NNLLO (next-to-ncxt-to leading logarithmic order) prediction in the 
Standard Model are needed. The virtual QED corrections at NNLO accuracy have been 
determined in a scries of articles quite recently [3-11]. A complete evaluation of the pho- 
tonic corrections covers additionally the real photon emission contributions and fermion pair 
production. 

In this talk, we discuss one class of Feynman diagrams for real photon emission, namely 
radiative loop corrections, 

e+ + e _ -> e+ + e" +7, (2) 

which are contributing at NNLO to reaction {T]). Their evaluation includes 5-point func- 
tions. Usually, the scalar, vector, and tensor functions of this type will be reduced to simpler 
one-loop functions. We also discuss an alternative approach, based on Mellin-Barnes repre- 
sentations of Feynman parameter integrals. 



*Presented by T.R. 

a A link to the slides of this contribution is [1] . See also [2] . 



LCWS/ILC2007 



2 Reduction of 5-point functions 



A critical point in an algebraic reduction of vector or higher tensor 5-point functions to scalar 

2- point, 3-point, and 4-point functions (in four dimensions) is the appearance of inverse 
Gram determinants. It is known that these inverse Gram determinants arc spurious [12-14] 
and that they may be canceled out in the final analytical expressions. For the approach 
proposed in [15], we have demonstrated this cancellation quite recently. Because that part 
of the presentation was decribed in some detail in other contexts [16,17], we don't repeat 
the material in these proceedings again. 

The focus will be on two questions: 

• Is the MB-approach useful for the numerical evaluation of the finite parts of scalar, 
vector, and tensor 5-point functions? 

• How to treat the infrared divergencies of these functions? 

3 Mellin-Barnes representation for massive 5-point functions 

The use of Mellin-Barnes (MB) integrals for the representation and evaluation of Feynman 
integrals has a long history, although a systematic use of it became possible quite recently. 
The replacement of massive propagators by MB-integrals was proposed in [18] for a finite 

3- point function. It was worked out for one-loop rt-point functions with arbitrary indices 
(powers of propagators) in d = 4 — 2e dimensions in [19-21], where also some of the related 
earlier literature is discussed, as well as the applicability to tensor integrals and to multi- 
loop problems. The aim was a replacement of massive by massless propagators. In [22], 
the Feynman parameter representation (or a-parameter representation, the difference plays 
no role here) was derived and then for the characteristic function of the diagrams an MB- 
represcntation was applied. Along this line, a systematic approach to MB-prcscntations for 
divergent multi-loop integrals was derived and solved for non-trivial massless and massive 
cases [23-29]. Since software packages like AMBRE.m [30] (in Mathematica, for the deriva- 
tion of MB-rcpresentations), MB.m [31] (in Mathematica, for their analytical and series 
expansion in e), and XSUMMER [32] (in FORM [33], for taking sums of their residues) 
became publicly available, quite involved integrals may be treated, see e.g. [34]. 

Of course, such a complicated task like the evaluation of - ideally - arbitrary Feynman 
integrals will not be finally solved with using one or the other method. In fact, already 
quite simple problems may be used to demonstrate the limitations of some approach. We 
will study here, with MB-integrals, some one-loop functions of massive QED as occurring 
in Bhabha scattering, with focussing on the 5-point function shown in Figure [1] 

Wc define 

with the chords Q;, 

di = (q- Qi) 2 - ml (4) 

This representation becomes unique after choosing one of the chords (and the direction of 
the loop momentum), e.g.: 

Q£ = 0, Q?=tf. (5) 
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The numerator A(q) contains the tensor structure, 

%) = {i,^,«Y,rtY,-}, (6) 

or may be used to define pinched diagrams; e.g. a shrinking of line 5 leads to a box diagram 
corresponding to 

h[ds] = e">* ( d \,,\, ■ (7) 
J nr d '' ! did 2 «3«4 

For details of the derivation of Fcynman parameter integrals we refer to any textbook on 
pcrturbative quantum field theory or to [35]. A Feynman parameter representation for Fig. 
HI is: 



I 5 [A(q)] = -e^ [ 1 l[dx j sfl-^xA ^±ll B(q), (8) 

■*° j=l \ i=l / * W 

with B(l) = 1, B(g") = Q»,B{q»q») = Q^Q" - \g^F{x)/{2 + e), and Q» = £ <BiQ?- The 
diagram depends on five kinematical invariants and the F-form in ([SJ is: 

F(x) = ml(x 2 + + X 5 ) 2 + [— s]xiX3 + [— V^S^Es + [— t]x2X4 + [— t']x2X5 + [— V2]xiX4,. 

(9) 

Henceforth, m e = 1. It is evident that the F-form cannot be made more compact. After 
the introduction of seven subsequent Mellin-Barnes representations, 



= -L [ dz[A(x)]'[Bx iXj ] 
J c 



_ R _ z T(R+z)T(-z) 



\A{x) + Bx iXj \ R 2m J c "" l " K *' J1 L ~*^ J T(R) (10) 

one for each additive term in F, we may perform the ^-integrations using a generalization 
of the integral representation of the Beta function: 

n T a a i~ l ^ r(ai) ■ • ■r(ajy) 
dxj x/ 5{l-xi x N ) = — -. 11 
■> r [a.\ H h a N ) 



LCWS/ILC2007 



The final MB-integral may be easily derived using our Mathematica package AMBRE.m 
[30] (see also examplel.nb and example2.nb of the package). The representation is five- 
dimensional after twice applying Barnes' first lemma. The integrals are well-defined on 
integration strips parallel to the imaginary axis, for a finite value of e = d/2 — 2. After an 
analytical continuation in e, preferrably done by MB.m [31], one gets a sequence of finite, 
multi-dimensional MB-integrals. We performed these steps and met the following situation 
for the terms proportional to 1/e and O(l): 

• scalar integrals: the MB-integrals are up to three-fold; 

• vector integrals: the MB-integrals are up to three-fold; 

• tensor integrals: the MB-integrals are up to five-fold. 

We performed some experimental calculations, but there is no need to go into more detail: 
A numerical evaluation of these integrals, especially of the five-dimensional ones, in the 
Minkowskian region, is not competative to the old-fashioned numerical packages like FF [36] , 
or LoopTools [36-38], which rely on the preceding algebraic reduction of all the 5-point 
functions to well-known scalar 2- to 4-point functions. 

For this reason, we restrict the discussion now to the infrared divergent parts only. As is 
well-known, they have a lower dimensionality, and here the MB-presentations are well-suited. 
As examples we will use the scalar and vector 5-point functions. 



4 Infrared singularities 



Let us consider first the scalar function. A set of five independent invariants may be read 
off from ©: 



s 
t 
t' 
V 2 
Va 



(Pi+Ps) 2 , 

(P4 +P5) 2 , 
(Pi +P2) 2 , 
2p 2 P3 ^^3, 
2p 4 p 3 ~ E 3 , 



and the two massless propagators are d§ = q 2 and d 2 = (q + p\ +Ps) 2 . In the IR-limit, 
where E3 — > 0, it will be t' fa t and < V 2 , V4 << s, \t\. The leading IR-singularities are 
easily found algebraically from the following decomposition: 



-1 



did 2 d 3 d4d 5 



2(g-Q 5 )(g-Q 2 ) l{2(q- Q 5 )(q - Q 3 ) 
V 2 



+ 



did 2 d3d4d5 
1 (2{q-Q 2 )(q-Q4) 



1 



1 



1 



did^dAds 
1 



did 2 d 3 d4 



d\d 2 ds d\d 3 d4 



d\d 3 d4 d\dAdn, 
(12) 



The 4-point functions depend on the variables (t,t',V 2 ) and (t, t', V4), respectively, and the 
leading IR-singularities of I5 trace back, by construction, to the two IR-divergent 3-point 
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functions: 



/ 



d d k i r d d k if d d k 1 



dxdidsd^d^ sVi J did^ds sVi J d\d 2 d 3 e 



F(t>) , F(t) 



sV 2 sV A 



(13) 



The integrals with numerators are constructed such that they are free of IR-singularities 
arising from the virtual photon lines. It is of importance here to observe that the denomi- 
nators V2 and V4 are proportional to the photon energy E3 and thus give rise to additional 
IR-problcms, stemming from the photon phase space integral over the squared sum of matrix 
elements; e.g.: 

d 3 p 3 A B(E 3 ) f u dE 3 _ 

VF~ F F / ~F~ - ln ^3j| - (14) 

Z£j 3 »3 Ji3 Jo ^3 

Here, one term (A/E3) comes from the real photon emission Born diagram, and the other 
one (B(Es)/ E3) from our pentagon diagram. After dimensional regularisation, this becomes 
evaluable and contributes also to the Laurent series in e. We learn from (fT4")l that a complete 
treatment of the IR-problem includes a careful control of the subleading (and in 4 dimensions 
non-integrable) terms like 1/Vi and \n(Vi)/Vj. This leads to phase space integrals with a 
behaviour like: 

dE 3 fa , . \ 2a + b c-2a\n{u) 



&ln(£ 3 ) + c) 



El~ d Ve v ' J 4e 2 2e 



cln(cj) + -{2a + b)\n{uj) + 0(e). (15) 



Evidently, one separates with the 3-point functions in (|13|) only a leading singularity, while 
we expect expressions like 



/ 



aict2«3«4«5 sV^e sVit sV 2 sVi sV 2 sv4 



Subleading singuarities may arise from the e-finite 4- and 3-point functions with pre-factors 
1/Vi. 

It is also evident that the whole above discussion immediately transfers over to vector 
and tensor integrals. 

Concentrating now on the IR-divergent parts, we may safely assume now the validity of 
the Born kinematics, including 

t' = t, (17) 

which is justified bcause of the vanishing photon momentum in this limit. This 'eats' another 
MB-integration (in the F-form @ one additive term vanishes), and the starting point of 
further discussions are four-dimensional MB-integrals. For the scalar pentagon: 



12 

+100+M; ~f r , 

if / dzii-sri-tri-v^ri-v,)- 3 -^-^ 3 -^ ^ 1 r 

{2m)' 1 il j r ri 3 ri4 

— ioc-\-Ui 

(18) 
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with a normalization Tq = T[—l — 2e], and the other T-functions are: 

ri = r[-*i], r 2 = r[-z 2 ], r 3 = r[-^ 3 ] J r 4 = r[i + z 3 ], 

r 5 = t[i + z 2 + z 3 ], r 6 = r[-* 4 ], r 7 = r[i + z 4 ], r 8 = r[-i - e - ^ - z 2 ], 

T 9 = r[— 2 - e - zi - Z2 — Z3 - za], Tio = F[—2 - e - Zi — Z3 - za], 
Tn = r[— e + 21 - z 2 + 24], Ti2 = r[3 + e + Zi + Z2 + Z3 + z 4 ], 

and, in the denominator: 

r i3 = T[—l — e — zi — Z2 — Z4], r i4 = T[-e - z\ — Z2 +za}. (19) 

The ^5 is finite if all T- functions in the numerator have positive real parts of the arguments; 
this may be fulfilled for finite e (here we follow the method invented in [25]): 

e = -| (20) 

The real shifts itj of the integration strips 7^ may be chosen to be: 

ux = -5/8, 

u 2 = -7/8, 

u 3 = -1/16, 

w 4 = —5/8, 

u 5 = -1/32. (21) 

The further discussion of the scalar case is very similar to that of the QED vertex function 
given in [39], so we may concentrate here on the results for the IR-divergent part: 



li R = l[ R (V 2 )+li R (V 4 ), (22) 
I s (V) 

H R m = ^^+/ s m)- (23) 



The explicit expressions for the inverse binomial sums solving the MB-integrals are obtained 
by applying the residue theorem (closing the integration contours to the left): 



2sV l e ^ 2n 

n=0 I 

n 



(2n + 1) 



with Iti being in accordance with (Q2|), and: 



h(V,Y = jvEt^-T^ [-21n(-V t )-3S 1 (n) + 2S,(2n + l)], (25) 



2 « (2„ + l, 



where we introduce the harmonic numbers Sk(n) = X)"=i l/* fc ' an d have to understand 
ln(-Fi) = \n(Vi/s) + ln[-(s + iS)/m 2 e }. 
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The series may be summed up in terms of polylogarithmic functions with the aid of 
Table 1 of Appendix D of [40]: 



jn - 

»=o ( 2n ) (2n + l) 



E 



o (i n )( 2 «+ 1 ) 



Si(n) 



E 



«=o ( 2n I (2n + l) 



5i(2n + l) 



(26) 



-4Li 2 (- 2 /)-41n( 2/ )ln(l + y) 



+ ln 2 (y)-2C 2 ] 



(27) 



v 2 -i 



2U 2 {y) - 4Li 2 (-I/) - 4In(») + I/) 



+21n(j/)ln(l-y) + -ln 2 ( ? /)-4C 2 



(28) 



with 



v = y{t) 



(29) 



For the vector and higher tensor 5-point functions one gets quite similar results. The 
IR-divergent pieces arise only from those contributions, which are proportional to the chords 
Q2 and Q5 of the massless internal lines (one of them is set to zero here, Q5 = 0): 



I V -l(V2,Vi) 



(30) 



The MB- integrals introduced in ^ will not get modified by the additional factors B(q^) 
etc., but the subsequent x- integrations will. For the vector integrals, we obtain: 



(31) 



and 



-\-ioo-\-Ui 



h(2) 



(2™) 4 Al 



-3 — €. — Zi—Z2 — Z3 — Z4 



-ioo~\-Ui 



12 

1 1 13 1 14 



(32) 



where it is F" = r ? with two exceptions: 



1 10 



K = rr-2ei 



r[-l - e - z\ - z 3 - zA 



(33) 
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After similar manipulations as described above, we obtain finally for the IR-divergent 
part of the vector pentagon (and, not discussed at all, the tensor pentagon): 

*6V] = Q% li R (V 4 ) + Q» U R (V 2 ) 7 (34) 

and 

H R kn = QZQ» 2 H R (V 4 ) + Q»Q" 5 4 R (V 2 ). (35) 

In the above derivations, we chose arbitrarily Q5 = 0. The leading and non-leading IR- 
divergent parts of the tensor functions are contained in those terms of the tensor decompo- 
sition, which are proportional to the chords of the massless internal lines, and they agree 
with the corresponding scalar functions. 

In conclusion, we have demonstrated, by analysing the loop functions without squaring 
matrix elements, that IR divergencies of scalar and tensor one-loop pentagon diagrams can 
be treated in a systematic, efficient way by using Mcllin-Barnes representations. The leading 
singularities of the e expansion of MB-integrals are obtained straightforwardly and have the 
same IR-structure as the vertex functions obtained by quenching. Both the leading and non- 
leading singular parts (the latter being kinematical end point singularities) can be expressed 
by a few well-known inverse binomial sums or, equivalently, polylogarithmic functions. The 
IR-structure of vector and tensor functions is completely reducible to that of the scalar 
function. 
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